The ^^N(p,7)^^0 reaction studied with a composite germanium detector 
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The rate of the carbon-nitrogen-oxygen (CNO) cycle of hydrogen burning is controlled by the 
^^N(p,7)^^0 reaction. The reaction proceeds by capture to the ground states and several excited 
states in ^^O. In order to obtain a reliable extrapolation of the excitation curve to astrophysical 
energy, fits in the R-matrix framework are needed. In an energy range that sensitively tests such 
fits, new cross section data are reported here for the four major transitions in the ^^N(p,7)^^0 reac- 
tion. The experiment has been performed at the Laboratory for Underground Nuclear Astrophysics 
(LUNA) 400 kV accelerator placed deep underground in the Gran Sasso facility in Italy. Using a 
composite germanium detector, summing corrections have been considerably reduced with respect 
to previous studies. The cross sections for capture to the ground state and to the 5181, 6172, and 
6792 keV excited states in ^^O have been determined at 359, 380, and 399 keV beam energy. In 
addition, the branching ratios for the decay of the 278 keV resonance have been remeasured. 

PACS numbers: 25.40.Ep, 25.40.Lw, 26.20.Cd, 26.65.+t 



I. INTRODUCTION 

The stellar rate of the carbon-nitrogen-oxygen (CNO) 
cycle of hydrogen burning [l|, Q is controlled by the 
slowest process, the i'^N(p,7)i^0 reaction [sj. In the 
Sun, hydrogen burning proceeds mainly by the compet- 
ing proton-proton chain, and the CNO cycle contributes 
only 0.8% to the energy production However, solar 
CNO hydrogen burning gives rise to neutrino emission 
lines from the /3+ decay of ^^N and ^^O It has re- 
cently been suggested [5] to use the expected CNO neu- 
trino flux data from the Borexino detector [6| and the 
planned SNO+ detector to measure the abundance of 
carbon and nitrogen in the solar core. This would ad- 
dress the so-called solar metallicity problem [sl, [9[ , which 
is given by the fact that the new solar metallicities [lo| 
lead to inconsistencies in the standard solar model. The 



correct interpretation of the expected CNO neutrino data 
requires, however, that the nuclear reaction rate of the 
CNO cycle, which is determined by the ^^N(p,7)^^0 cross 
section, be known with sufficient precision. 

The ^^N(p,7)^^0 cross section cf{E) can be parameter- 
ized using the astrophysical S-factor 



S{E) = aEexp \212A/Ve 



(1) 
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with E denoting the energy in the center of mass system 
in keV. 

The excitation function has been studied previously 
pj]-[l3[ e.g.], and these data determine the recommended 
value in the current nuclear reaction rate compilations for 
astrophysics p^-[T6| . Subsequently, a number of new ex- 
perimental and theoretical results on this reaction have 
been reported p!7l-[26|, showing that the recommended 
value of the reaction rate [ij-llq has to be revised down- 
ward by a factor of two. In particular, capture to the 
ground state in ^^O (fig. [1]) was shown to be strongly 
suppressed (T^-Hsl. This reduction is now adopted in a 
very recent compilation p7|. However, some open ques- 
tions remain. 

In particular, two groups have in recent years pre- 
sented cross section data and performed R-matrix fits 
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FIG. 1: Level scheme of ^^O, in keV [H, [2^. The most 
important 7-transitions are denoted by arrows. 

based on their new data: LUNA [111, ^ and TUNL 
[22i] . These two works show excellent agreement when 
it comes to the most important contribution to the to- 
tal S-factor, namely capture to the state at 6792 keV: 
1.20±0.05 keVbarn [23] and 1.15±0.05 keVbarn [H, 
respectively. However, their results differ by much 
more than the quoted uncertainties when it comes to 
the second most important contribution, capture to 
the ground state: Whereas LUNA reported^ *5'gs(0) = 
0.25±0.06 keVbarn jpl, the TUNL value is double that, 
0.49±0.08 keVbarn [22|. This discrepancy amounts to 
about 15% of the total extrapolated S'tot(O), dominating 
the uncertainty. 

The only significant methodical difference between the 
two fits from LUNA [2l|, Hsj and TUNL [22| is the treat- 
ment of high-energy data. The LUNA fit is a global fit, 
based on the LUNA data presented in the same paper, 
and on the Schroder et ai data [l3| which had been 
corrected for the summing-in effect. The TUNL fit, on 
the other hand, is a partial fit based solely on the TUNL 
data presented in the same paper, with the higher-energy 
R-matrix poles kept fixed based on a previous fit of the 
Schroder et ai data . The starting values and general 
procedure for both fits are otherwise the sa me il8i]. 

The experimental data points by LUNA [211, |23| and 
TUNL [22I are generally in agreement with each other, 
but they show some systematic uncertainty due to the 
fact that both groups had employed large germanium 
detectors in close geometry. This arrangement had been 
chosen in order to obtain a high enough detection effi- 
ciency for the weak ground state capture line. However, 



in this way both groups also incurred true coincidence 
summing-in corrections of more than 100% for the ground 
state data. Such a large correction entails considerable 
systematic uncertainty. 

The aim of the present work is to address the conflict- 
ing extrapolations [2l|, [22| in two ways. The experimen- 
tal problem of the previous high summing-in correction is 
solved by using a Clover detector. The problem of the se- 
lection of the database is solved by providing the ground 
state cross section relative to that for the well-known cap- 
ture to the state at 6792 keV. The present relative data 
can then be added to one particular data set without in- 
troducing additional scaling uncertainty. Alternatively, 
they can be rescaled to absolute data using an overall fit 
of 6792 keV capture based on several independent works, 
strongly reducing the scaling uncertainty. 

For the present experiment, the energy range of E 
= 317-353 keV has been selected, far enough above the 
259 keV resonance to limit resonant contributions, and at 
the same time a region where a sensitive minimum [l8| 
of R-matrix fits is observed. In principle, such a mea- 
surement would also have been possible at £^ ~ 170keV, 
in a second sensitive minimum. However, the yield is a 
factor 100 lower there, so that the present energy range 
was chosen for practical purposes. 

The present relative cross section data have been pub- 
lished previously in abbreviated form The present 
work provides full details of that experiment and analy- 
sis. In addition, new branching ratios for the decay of the 
259 keV resonance obtained in even farther geometry are 
presented here. The absolute off-resonance ^^N(p,7)^^0 
cross section for capture to the ground state and the 5181, 
6172, and 6792 keV excited states is derived at Ep = 
359, 380, and 399 keV. In order to improve the reliabil- 
ity, this latter analysis is performed in two independent 
ways, namely by the 7-line shape method [29| and by the 
classical peak integral approach. 



II. EXPERIMENT 

The experiment was performed in the Laboratory for 
Underground Nuclear Astrophysics (LUNA) at the Gran 
Sasso National Laboratory (Italy) [30|,[3l|. At the LUNA 
site, the 7-ray laboratory background for E^ > 3 MeV is 
strongly reduced due to the rock overburden equivalent 
to 3800 meters water [32,^]. Also for E^ < 3 MeV with 
proper shielding the 7-ray background has been found 
to be much lower than in comparable laboratories at the 
surface of the Earth [s^ • The unique location of LUNA 
has enabled the study of several nuclear reactions of as- 
trophysical importance [2l|, [35| - [4Q| . 



A. Setup 

^ Si(0) denotes the S-factor, extrapolated to zero energy, for cap- 1 — 1 

ture to the state at i keV in ^^O. 5gs(0) and ^tot(O) refer to The LUNA2 400 kV accelerator [29^ provided a H+ 

ground state capture and to the total S-factor, respectively. beam of Ep = 359, 380, and 399 keV, with 0.25-0.45 mA 
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FIG. 2: (Color online) Schematic view of the experimental 
setup. 
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FIG. 3: (Color online) Scans of the Ep = 278 keV resonance 
at the beginning of the experiment (black squares), one day 
with 29 C dose later (red triangles) , and at the end of the 
experiment after a total accumulated dose of 267 C (blue cir- 
cles) . 



intensity. The ion beam passed a collimator of 5 mm di- 
ameter, which absorbed a few percent of the full beam in- 
tensity, and a cold trap cooled by liquid nitrogen (fig. [2]), 
before hitting the target. Secondary electrons emitted 
from the target surface were suppressed by applying - 
300 V suppression voltage to the cold trap. The repro- 
ducibility of the current from run to run is estimated to 
be 2%. 
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FIG. 4: (Color online) (a) 7-ray detection efficiency for the de- 
tector at 9.5 cm distance from the target, as determined with 
radioactive sources and the two-line method at the 259 keV 
resonance. Solid (dashed) curve, efficiency for addback (sin- 
gles) mode, (b) Residuals. The data point at 7556 keV is not 
yet corrected for summing-in. It was excluded from the fit, 
and is shown here for illustration only. The pair of 7-rays 
at 5181 and 2375 keV was also not included in the fit but is 
plotted here as a check on the reliability of the curve. 



under 55° angle. In order to obtain its stoichiometry, 
the stopping power at the resonance energy, the beam 
current and the strength of the monitor resonance must 
be known. For the stopping of protons in titanium and 
nitrogen, the values from the SRIM software [4l| have 
been used. For the strength of the resonance, ujj = 
13.1±0.6meV was adopted, the recommended value from 
Ref. [23. Based on this number, a stoichiometric ratio 
Ti:N of 1:0.93 has been determined. The target stoi- 
chiometry gives rise to 6% systematic uncertainty in the 
absolute cross section results, mainly from the reference 
ujj value. 

In order to properly correct for the change of the tar- 
get under intense proton bombardment, during the ex- 
periment the target profile was monitored every day by 
scanning the = 278 keV resonance (fig. [3]). The sharp 
low-energy edge of the profile is given by the convolu- 
tion of the O.lkeV energy spread of the beam [29] and 
the 1.06 keV natural width of the resonance [42]. On the 
ensuing constant plateau, the step height is proportional 
to the inverse of the effective stopping power per ^^N 
nucleus in the compound. 



B. Target 

A titanium nitride target produced by reactive sput- 
tering at Laboratori Nazionali di Legnaro was used for 
the experiments. The target had 60keV energetic width 
at the Ep = 278 keV resonance (fig. [3]), when irradiated 



A reduction of up to 7% in the integral of the target 
profile was observed from day to day, with a typical pro- 
ton dose of 24 C (1.5 • 10^° H+ ions) deposited on the 
target per day. It is estimated that the target composi- 
tion is known with 5% precision for any given time during 
the experiment. 
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C. Detection of emitted 7-rays 

The 7-rays emitted from the target were detected in a 
Eurisys Clover-BGO detection system [43| placed at an 
angle of 55° with respect to the beam axis. The front end 
of the Clover detectors was at 9.5 cm distance from the 
target. For the branching-ratio measurement (sec. lIII D]) . 
the front end was placed at 19.5 cm distance from the 
target instead. 

The output signal from each of the four Clover crystals 
was split into two branches called branch 'S' and branch 
'A'. For branch 'S', each of the four signals was amplified 
and digitized separately, and the four spectra were gain- 
matched and summed in the offline analysis, giving the 
so-called singles mode. 

For branch 'A', the preamplifier output signals were 
gain-matched and added in a homemade analog summing 
unit. The added signal was then amplified and digitized, 
giving the so-called addback mode spectra. Typical res- 
olutions for addback (singles) mode were 9keV (3.3 keV) 
at 1.3 MeV and 12keV (6keV) at 6.8 MeV. For experi- 
ments off the 259 keV resonance, the addback mode data 
were recorded in anticoincidence with the BGO escape- 
suppression shield to reduce the Compton background. 

The 7-ray detection efficiency was measured using 
-^^^Cs and ^^Co radioactive sources calibrated to 1.5% 
and 0.75% (la confidence range), respectively. The ef- 
ficiency curve (fig. 21 upper panel) was then extended 
to high energy based on spectra recorded at the 259 keV 
1/2+ resonance, using the known 1:1 7-ray cascades for 
the excited states at 6172 and 6792 keV [42]. The 7- 
ray s from the decay of this 1/2+ resonance are isotropic 
[42|. The angular correlations of 8-10% between priniary 
and secondary 7-ray are experimentally well known [4^ . 
They result in up to 0.4% correction on the efficiency 
curve, because they affect the summing-out correction. 
For the worst case, the 6792 keV 7-ray, the calculated 
summing-out correction is 3.6% in addback mode (1.1% 
in singles mode), with an assumed relative uncertainty 
of 20%. This result is consistent with a GEANT4 [ig 
simulation showing (4.5±1.8)% correction. 

As a check on the quality of the efficiency curve, the 
experimental cascade ratio for the 5181 keV excited state 
(not used in the fit) was found to be reproduced within 
1% statistics (fig. [H lower panel), again assuming 1:1 7- 
ray cascade ratio [42]. 



III. DATA ANALYSIS AND RESULTS 

The first part of the analysis concentrated on the ratio 
of the cross sections for radiative proton capture to the 
ground state and the fourth excited state at 6792 keV 
in ^^O, determined with the detector at 9.5cm distance 
from the target. These relative data have been reported 
previously in abbreviated form [25| and are discussed in 
details in sec. IIII A[ Subsequently, also absolute cross 
section data for the four most important 7-transitions 



are derived from the spectra. This analysis is performed 
both by classical peak integrals for the addback mode 
data (sec. IIII Bp . and by 7-line shape analysis for the 
singles mode data (sec. IIII Cp . Finally, by moving the 
detector to 19.5 cm distance from the target, more precise 
branching ratios for the decay of the 259 keV resonance 
are presented (sec. IIIID]) . 



A. Ratio of the cross sections for capture to the 
ground state and the 6792 keV excited state in ^^O 

For the relative analysis, the number of counts in the 
ground state capture peak at ^ 7600 keV is com- 
pared with the number of counts in the secondary 7-ray 
at 6792 keV (fig. [5|). In such an analysis, only the rela- 
tive uncertainty when extending the efficiency curve over 
this limited energy range contributes to the uncertainty 
of the ratio (0.8% effect). 

The 6792 keV counting rate contains some on-resonant 
contribution. This is due to the 60 keV (full width at half 
maximum) thick target. When the beam slows down to 
the strong resonance at Ep = 278 keV, it still finds some 
TiN in the tail of the target. In order to correct for 
this effect, the primary 7-rays for capture to this level 
are analyzed, as well, and the 6792 keV counting rate 
is rescaled with the resonant /off-resonant ratio as ob- 
tained from the low-energy primaries (fig. [6]). The reduc- 
tion in 6792 keV counting rate by the escape-suppression 
shield contributes 1.2% to the final uncertainty, and the 
summing-out correction for this peak contributes 0.6%. 

Based on these data, the ratio 



^GS/6792(^) 



^GS(^) 
(J6792{E) 



(2) 



has been calculated (table H]). The present data supersede 
the data published previously in abbreviated form [25j . 
due to an upgraded background determination (fig. [6l 
blue dashed lines), described in section UlIBI The ratio 
depends only on the counting rates for the E^ ^ 7600 keV 
ground state capture 7-ray, for the E^ = 6792 keV 7- 
ray (corrected for resonant capture as described above), 
and on the ratio of the 7-detection efficiencies at E^ ^ 
7600 and 6792 keV. For the ground state capture 7-ray, a 
summing-in correction of up to 30% (4.3%) for addback 
(singles) mode was taken into account (table HI last col- 
umn). 

When computing ^^08/6792(^)7 the current measure- 
ment and the target stoichiometry and profile cancel out, 
eliminating the major sources of uncertainty. There- 
fore, the relative analysis method allows to derive data 
with much better precision than for absolute data. The 
present relative data can then be rescaled with averaged 
data for the well-studied cross section for capture to the 
6792 keV state, and uniquely precise data for capture to 
the ground state can be obtained. 

The effective interaction energies have been deter- 
mined for each 7-line with two methods: First, the cen- 
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FIG. 5: (Color online) High-energy part of the 7-ray spectra for addback mode (red full line) and singles mode (green dashed 
line), (a)-(c): Off-resonant spectra at Ep — 399, 380, and 359 keV, respectively, with the detector at 9.5cm distance from the 
target, (d): On-resonance spectrum on the = 278 keV resonance, with the detector at 19.5 cm distance from the target. 
The resonant contribution by the tail of the target is well visible also in the off-resonant spectra. The contaminant peaks stem 
from the "^^F(p,a7)"^^0 reaction. 



TABLE I: Cross section ratio i?GS/6792(^) and relative un- 
certainty. The size of the summing-in correction is also given. 
The present data supersede Ref. [2^, due to an improved 
background determination. 



E [keV] 


mode 


^GS/6792(^) 

[10-2] 


stat. syst. 


Summing-in 






uncertainty 


correction 


315.3±1.3 


addback 


5.24 


11% 5.4% 


30% 




singles 


5.22 


15% 2.7% 


4.3% 


333.1±1.0 


addback 


5.33 


4.8% 3.9% 


21% 




singles 


5.58 


11% 2.5% 


3.4% 


353.3=bl.O 


addback 


5.20 


3.5% 3.5% 


19% 




singles 


5.43 


8.0% 2.3% 


3.2% 



troid of the off-resonant primary 7-line has been used, 
taking into account the reaction Q-value and 7-level en- 
ergies. Second, the average energy, weighted with the 
predicted counts from the known target profile and the 
expected energy dependence of the cross section from the 



R-matrix S-factor curve The two values were never 
more than 2.6 keV apart, and their average was adopted 
for each line. The results are slightly different for ground 
state capture and capture to the 6792 keV state, because 
the S-factor curve from previous R-matrix fits has a dif- 
ferent slope for these two transitions. Therefore, the av- 
erage of the two values is adopted as effective energy to 
be connected with the cross section ratio ^^08/6792(^)7 
with the assigned la error bar covering both effective 
energy values. 

For the relative data, the total systematic uncertainty 
is 3.5 - 5.4% in addback mode (table |TT1). For singles 
mode, due to the lower summing corrections, it is 2.3 - 
2.7%. 

For all three data points, the addback and singles 
mode data are in good agreement. Due to the higher 
7-efficiency of the addback mode data (which, in turn, is 
due to the well-known addback factor of Clover-type de- 
tectors [isl , which has been redetermined for the present 
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detector and geometry [33|) and due to the background 
reduction achieved by the escape-suppression shield for 
the addback mode data, the addback data have much 
better statistics than the singles mode. Therefore, the 
addback data are adopted for the further analysis despite 
their slightly higher systematic uncertainty. 



B. Absolute cross sections based on the peak 
integrals of the addback mode data 

As a second step, the absolute cross section for capture 
to the excited states at 5181, 6172, and 6792 keV and to 
the ground state of ^^O has been derived, accepting that 
the systematic uncertainty (table (III) includes now also 
the contributions from current measurement, target sto- 
ichiometry and profile, and absolute detection efficiency. 
Only the addback mode data were considered. 

In order to obtain the net counting rate, a straight-line 
background based on two flat regions to the left and right 
of the region of interest (ROI) has been subtracted from 
the integral over the ROT This procedure was applied 
for every secondary except for the decay of the 6172 keV 
excited state, where a different method was applied. It 
was repeated for each transition of the run at 399 keV, 
both for the primary (resonant and non-resonant) and 
secondary 7-rays. 

However, in many cases it was not possible to apply 
this method of background determination: At = 359 
and 380 keV, the off-resonant part of the primaries lie 
close to the resonant peak (fig.O second and third row). 
The secondary at ^ 6172keV was problematic, as 
well, due to the ^^F(p,a7)^^0 background peak at E^ ^ 
6130 keV (fig. [SJ. For these spectra, a different method 
was instead used to estimate the background: The ra- 
tio between the difference in average counts per channel 
observed to the left and right of the peak, and the net 
area of the peak itself, was calculated. The ratios ob- 
served on the resonance, where no additional resonant 
contribution exists and where beam-induced background 
is negligible, have then been used to calculate the back- 
ground at the same 7-energy in the problematic spectra. 
For those problematic spectra, a minimum uncertainty of 
5% has been assumed for the quantity subtracted from 
the raw integral of the ROI. Finally, it was ensured that 
the 1(7 uncertainty of the counts includes also results with 
different choices of background regions. 

The net counting rate was then determined from the 
secondary 7-ray, rescaled for its non-resonant /resonant 
contributions determined by the primary 7-rays. Based 
on the counting rate, the target stoichiometry and profile 
(sec. Ill Bl fig. [3]), the beam current measurement, and 
the 7-detection efficiency (fig. |4|), the cross section was 
calculated for these transitions. The angular distribution 
was assumed to exhibit negligible contributions from all 
Legendre polynomials except for zero and second order. 
The second order Legendre polynomial cancels out at the 
present detection angle of 55°. 



TABLE II: Systematic uncertainties affecting cross section ra- 
tios ("relative", sec. IIII A)) and absolute cross sections ("ab- 
solute", sees, nil 51 and UlTC]) for addback mode. 



Affecting data... Description Amount 

Relative Absolute Summing-in, ground state line 3-5% 

Relative Absolute Escape-suppression efficiency 1.2% 

Relative Absolute Slope of 7-efiiciency curve 0.8% 

Relative Absolute Summing-out 0.6% 

Absolute Target, original stoichiometry 6% 

Absolute Target, profile change 5% 

Absolute Assumption on S-factor slope 1-9% 

Absolute Beam current reproducibility 2% 

Absolute Normalization of 7-efficiency 1.8% 

Relative Total, addback mode 3.5-5.4% 

Absolute Total, addback mode 9-12% 



For the determination of the astrophysical S-factor 
from a single data point, it is necessary to make some 
assumption on the relative shape of the S-factor curve. 
For the present analysis, the S-factor was assumed to vary 
over the target thickness as given by the previous LUNA 
R-matrix curve [23|. In order to check the uncertainty 
introduced by this assumption, the present analysis was 
repeated assuming a flat S-factor, and the full difference 
(1-9%, depending on the transition and beam energy) 
was adopted as s yste matic uncertainty. The effective in- 
teraction energy [43 was calculated based on the known 
target profile and the assumed S-factor behaviour. The 
uncertainties are half of the difference obtained by using 
a flat S-factor instead of the LUNA's curve (23| . 

C. Absolute cross sections based on the 7-line 
shape analysis of the singles mode data 

Subsequently, the absolute cross section for capture to 
the excited states at 5181, 6172, and 6792 keV and to 
the ground state of ^^O has also been calculated based 
on the 7-line shape analysis approach. To this end, only 
the singles mode data, which are essentially free from 
summing corrections, have been used. This approach is 
thus complementary to the one described in the previous 
section, which calculated peak integrals and used only 
the addback mode data. 

The 7-line shape analysis method has been described 
previously in details [21, 23, 29], so it will only be out- 
lined here. The analysis of the line shape of the primary 
7-ray is possible because the observed line shape of a 
primary transition is determined by the cross section be- 
havior (t{E) in the proton energy interval spanned by the 
incident beam during the slowing-down process in the 
target. Each center-of-mass beam energy E (at which 
the reaction takes place) corresponds to a 7-ray energy 

E^ = E -\- Q — E^ -\- A£^Doppler — A£^Recoil, (3) 

with Q the reaction Q- value, E^ the energy of the excited 
state, and A£^Doppier/Recoii the appropriate Doppler and 
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FIG. 6: (Color online) Low-energy part of the 7-ray spectra for addback mode (red full line), after subtraction of the laboratory 
background. For completeness, the laboratory background is also included (black dotted line). Rows from top to bottom: 
= 399, 380, 359 keV. Columns from left to right: Primary 7-ray for capture to the excited state at 6792, 6172, 5181 keV. The 
peak from resonant capture (by the tail of the target) is clearly visible at the left of each panel. The non-resonant capture 
has a shape reflecting the profile of the target, convoluted with the energy-dependent cross section. The analysis of the 6172 
primary (central column) is hampered by the strong laboratory background line; this is reflected in higher uncertainty for 
this transition. The regions of interest and the assumed background are shown by blue dashed lines. 



recoil corrections. The 7-line shape is also influenced by 
the energy loss of the protons in the target, because the 
stopping power of the protons in titanium nitride is a 
function of proton energy [41]. 

The number of counts in channel i of the 
7-spectrum, corresponding to the energy bin 
[E^i^E^i-\-SE^]^ where SE is the dispersion in units 
of keV per channel, is given by the expression 

_ cT{Ei)6E^r]fe{E^i)bk 
e{E,) 

for Ei < E. Here Ei is the center-of-mass proton en- 
ergy corresponding to channel E is the incident proton 
energy in the center-of-mass, cr{Ei) is the cross section 
under study, r]fe{Ej^i) is the 7-ray detection efficiency, 
£{Ei) is the stopping power and b/e is the branching of 
the transition under study. The conversion from ^ 
to includes the Doppler and recoil effects, as shown in 
eq. (|3]) . The resulting count rate is folded with the known 



energy resolution of the 7-ray detector to obtain the 
experimental line-shape. 

To facilitate the fit, the cross section (t{E) entering 
into eq. (j4]) is then parameterized, in the limited energy 
window defined by the target thickness A^^xarget, as the 
sum of a resonant term described by the Breit-Wigner 
formula, and a non-resonant term, for which a constant 
astrophysical S-factor 5'nr is assumed: 

Here, A is the de Broglie wavelength, ujj the strength 
value of the 259 keV resonance (here, 12.9 meV was 
used [23I, very close to the recently recommended value 
of 13.1 me V [23), the energy of the resonance, F 
the energy-dependent total width of the resonance, and 
g-27rr7 Sommcrfcld parameter. Since the branching 

ratios and the ujj of the resonance are kept fixed, the 
free parameters in this procedure are the non-resonant 
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FIG. 7: (Color online) Typical 7-ray line shape obtained at = 380 keV, in singles mode, for capture to the 6792 keV state. 
The dashed black line is the assumed background, fitted outside the peak area. The dash-dotted (green) line corresponds to 
the expected resonant contribution (first part of eq. (O), and the dotted (blue) line to the fitted non resonant part (second part 
of eq. (O). The solid (red) line is the sum of these two components and the background. 



S-factor 5'nr, the background parameters and the energy 
of the beam. They are fitted to best reflect the shape of 
the primary 7-line by reducing the x^- After the fit has 
converged, the cross section a under study here is given 
by the average of the a{Ei) values, weighted for their 
contribution to the total statistics. 

Figure [71 shows a typical case for the primary 7-ray 
spectrum, together with the fit described above. The 
drop in the 7-ray yield towards lower energies reflects 
mainly the drop of the cross section due to the lower 
Coulomb barrier penetrability at lower energy. The en- 
ergy of the high energy edge of the peak provides an inde- 
pendent cross-check on the assumed beam energy from 
the accelerator energy calibration [29|. Possible varia- 
tions of the stoichiometry of the titanium nitride target 
during the beam bombardment have been monitored as 
described above fsec. lIIBl and fig. [3]). 

The final astrophysical S-factor obtained from the line- 
shape analysis described in the present section was found 
to be in excellent agreement with the data from the peak- 
integral approach described in the previous section. It 
should be noted that while the present line-shape analy- 
sis is based on the singles mode spectra, the peak inte- 
gral analysis is based on the addback mode data. The 
agreement between these two approaches confirms their 
reliability. 

The final S-factor values from the present experiment 
are obtained by forming the simple average value of the 
two approaches (secs. lIII Bl and HIl Cp . The data are sum- 
marized in tableHni and plotted in fig.[8l 



D. Branching ratios for the decay of the 259 keV 
resonance, obtained in far distance 

In order to determine the branching ratios for the de- 
cay of the 259 keV resonance {E^ = 7556 keV in ^^O), 
the Clover detector was moved to a farther geometry, 
with its front face at 19.5 cm distance from the target 
position, again at an angle of 55° with respect to the 
beam direction. For the branching ratio analysis, both 
addback and singles mode data have been analyzed and 
were found to agree within their statistical uncertainty 
in all cases. In the following text, only the singles mode 
data will be discussed. 

The detection efficiency was again established as de- 
scribed above (sec. Ill C]) . with an analogous quality of the 
efficiency curve as the one shown for the 9.5 cm geometry 
(fig-Hj). It should be noted that the efficiency curve does 
not depend on the branching ratios, just on the assump- 
tion of 1:1 cascade ratios without feeding or intermedi- 
ate decay corrections for the two transitions through the 
states at 6172 and 6792 keV, and on the assumption of 
isotropy 

For the determination of the decay branchings of the 
259 keV resonance, only the secondary 7-rays at 5181, 
5241, 6172, and 6792 keV and the ground state pri- 
mary 7-ray at 7556 keV were used (fig. El bottom panel). 
Therefore only the relative 7-efficiency in the limited en- 
ergy range 5181-7556 keV is needed. Owing to the good 
quality of the 7-efficiency curve, over this limited energy 
range the efficiencies relative to the 6172 keV normaliza- 



9 



0.01 




200 
E[keV] 



200 
E[keV] 



FIG. 8: (Color online) S-factor for capture to the excited states at 5181, 6172, and 6792 keV and to the ground state in ^^O. 
Data: Black diamonds [l3], blue squares [2ll . [23|. green circles [S^], red full triangles (present work, average of sees. IIII Bl and 
IIII C|) . R-matrix fits: Black dotted curve |l9l ]. blue dash-dotted curve [2ll. [23|. green dashed curve [23], black full curve [27], red 
long-dashed curve [2^. — For ground state capture, the black inverted triangles represent the present relative data f sec. IIII A)) , 
rescaled with the averaged S-factor for capture to the 6792 keV state as described in the text. For capture to the 5181 keV 
state, no R-matrix fits are given in Refs. p^, [2^. Error bars reflect the statistical uncertainty. 



TABLE III: S-factor results for capture to the ground state and to the excited states at 5181, 6172, and 6792 keV. The effective 
energy E is given in keV, the S-factor S in keV barn, and the relative uncertainties for S in percent. 



Capture to ground state 

E SgS Astat Asyst 


Capture to 5181 keV state 

E S'siSl Astat Asyst 


Capture to 6172 keV state 

E aS'6172 Astat Asyst 


Capture to 6792 keV state 

E *S'6792 Astat Asyst 


314.6±1.0 0.074 11% 12% 
333.6±1.0 0.061 5% 11% 
353.9±1.0 0.061 4% 10% 


310.6±2.2 0.370 16% 11% 
327.6±1.6 0.218 12% 12% 
350.9±2.5 0.128 13% 10% 


310.5±1.0 1.072 8% 12% 
326.6±1.0 0.406 18% 12% 
351.1±2.2 0.220 15% 10% 


315.9±1.3 1.495 5.0% 9% 
332.6±1.0 1.245 3.0% 9% 
352.7±1.0 1.157 1.7% 9% 



tion point are known on the level of ±0.5%, enabling a 
precise determination of the branching ratios. 

For the major transitions through the excited states at 
5181, 6172, and 6792 keV, the present branching ratios 
(tab. lIVp are in excellent agreement with the modern 
literature [22|, [23|. However, some minor discrepancies 
arise when it comes to the minor transitions. 

The ground state transition has been the subject of 
discussion in recent years. It is now well-known that 
the previously accepted value of (3.5±0.5)% [III, EH, Esj 
was much too high, probably due to summing-in. The 



two most recent previous branching ratio measurements 
[22I . [23! were both performed at about 20 cm distance, 
where there is still more than 10% summing-in correc- 
tion. The present value of (1.49±0.04)% has been ob- 
tained at 19.5 cm distance, with just 2.0% summing-in 
correction for the singles mode data, much less than 
in previous works. Note that the value (1.53±0.06)% 
from an abbreviated version of the present work [25] had 
been obtained in closer geometry, at 9.5 cm distance, with 
7.4% summing-in correction. The present (1.49±0.04)% 
ground state branching supersedes all previous LUNA 
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TABLE IV: Branching ratios for the decay of the 259 keV 
resonance {E^ = 7556 keV in ^^O) obtained with the Clover 
detector in singles mode, at 19.5 cm distance from the target. 
The numbers are compared with previous data [2^ . [23l . . 

Branching [%] 
Ajzenberg- TUNL |22J LUNA [23] LUNA, 
Selove [42] present work 



7556^ 


3.5=b0.5 


1.70±0.07 


l.GdzO.l 


1.49di0.04 


^5181 


15.8±0.6 


17.3±0.2 


17.1±0.2 


17.3±0.2 


^5241 






0.6db0.3 


0.15db0.03 


^6172 


57.5±0.4 


58.3±0.5 


57.8=b0.3 


58.3±0.4 


^6792 


23.2±0.6 


22.7±0.3 


22.9=b0.3 


22.6±0.3 



0.20=b0.05 keV barn. That value is lower than the recently 
recommended 0.27±0.05keVbarn i27|], but still in agree- 
ment given the error bars. The difference is mainly due 
to the fact that in the present work, only the present and 
the Schroder ^] data (corrected for summing-in ^]) are 
included. The data from Refs. [22|, [23| are excluded due 
to concerns about the summing corrections. In Ref. p7| , 
instead, the data from Refs. |g2|, i23| have also been in- 
cluded in the fit. 



V. SUMMARY AND OUTLOOK 



branching ratio measurements of the 259 keV resonance, 
i.e. [23, 25]. 

For the transition to the 5241 keV state, the previous 
(0.6±0.3)% value [23| was possibly affected by feeding 
through higher- lying excited states. Based on the differ- 
ence between 5241 and 7556 5241 7-rays, this 
feeding contribution amounts to (0.20±0.10)% of the to- 
tal decay branching. It is probably due to the 6859 keV 
state, which decays to 100% to the 5241 keV state (42j . 
However, such a weak feeding could possibly also arise 
through the 6172 or 6792 keV states, so in absence of 
conclusive evidence this (0.20±0.10)% is not assigned to 
any transition. 

For the transition to the 5181 keV state, the present 
data confirms the slightly higher modern values |22l , [23| 
with respect to the compilation [42|. 



IV. R-MATRIX ANALYSIS OF GROUND 
STATE CAPTURE 

For the purpose of an R-matrix analysis, the present 
relative data (sec. IIII A( table H]) have been renormal- 
ized using a weighted average S-factor for capture to the 
6792 keV state. Based on these values and the corrected 
Schroder data [l3|, [2l|, a new R-matrix fit for ground 
state capture has already been presented in the abbrevi- 
ated form of the present work [25| . The present updated 
relative data are close to the values published in abbrevi- 
ated form [25I , so this update does not warrant a revised 
fit. 

Also the present absolute data (sec. IIII Bp do not sig- 
nificantly deviate from the relative data, renormalized as 
stated above (fig. [51 bottom right panel). It should be 
noted that the present absolute data for capture to the 
6792 keV state (fig. [HI bottom left panel) are in excellent 
agreement with previous data (22I. |23| and R-matrix fits 
(22I . [23! , confirming that the renormalization procedure 
was adequate. By design the absolute data have higher 
uncertainty than the relative data (table [TT]) that have 
already been included in the fit |25| , so no new R-matrix 
fit is attempted here. 

The previous fit [25j is instead shown again here 
(fig. [HI bottom right panel), leading to S'gs(O) = 



The ^^N(p,7)^^0 reaction has been studied with a 
composite Clover-type detector at the LUNA under- 
ground facility at = 359, 380, and 399 keV, in an en- 
ergy range important for future R-matrix fits of capture 
to the ground state in ^^O. Precise cross section ratios for 
ground state capture relative to capture to the 6792 keV 
state have been presented, updating and extending their 
previous abbreviated publication |25| . 

The present, precise relative cross section data (table [I]) 
helped resolve the discrepancy between the previo us, con- 
flicting extrapolations for ground state capture [2l|, [22| , 
in favour of Ref. [2l| . The present recommended value of 
5'gs(0) = 0.20±0.05keVbarn is based on a dataset where 
the summing-in correction is not larger than 50% [l3| for 
the high-energy data and not larger than 30% for the 
present, lower-energy data. 

The present absolute cross sections for capture to 
the excited states at 5181, 6172, and 6792 keV (fig. [HI 
table imp have been obtained with two independent 
analysing methods (sees. IIII Bl and [TTTC]) . They are gen- 
erally in good agreement with previous works (22I, |23| 
and in some cases more precise. They are in overall good 
agreement with the most recent R-matrix fit [l^- Be- 
cause of their limited energy span, the present data alone 
cannot form the basis of new extrapolations. However, 
they may serve as useful reference points in an energy 
range that may be accessible not only at LUNA, but also 
at future underground accelerators. 

The new branching ratios for the decay of the 259 keV 
resonance that are shown here improve the precision of 
the database for this resonance. Since this resonance is 
often used as normalization point for experimental work 
on the ^'^N(p,7)i^0 reaction \2^^, e.g.], this im- 
proved information again facilitates future precision stud- 
ies of this reaction. 

The present data are an important ingredient in up- 
dates of the standard solar model [5i, ^8^, iSj. When ex- 
perimental data for the flux of solar CNO neutrinos due 
to the /3-decay of ^^N and ^^O become available from 
Borexino ^ or SNO+ 0, precise ^^N(p,7)-'^^0 cross sec- 
tions may contribute to a direct measurement of the solar 
metallicity through a comparison of CNO and ^B neu- 
trino fluxes [5[. 

Possible next steps in improving the precision for the 
extrapolated S-factor of this reaction [23 are to re-study 
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the cross section at higher energies [265, order to im- 
prove the extrapolation, and a remeasurement of the 
strength of the 259 keV resonance. 
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